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ERIC LUgON^ 

Abstract. We investigate the long-time asymptotics of the fluctuation SPDE in the Ku- 
ramoto synchronization model. We mainly establish the linear behavior for large time and 
weak disorder of the quenched limit fluctuations of the empirical measure of the oscillators 
around its McKean-Vlasov limit. This is carried out through a precise spectral analysis of 
the unbounded evolution operator, using arguments of perturbations of self-adjoint operators 
and analytic semigroups. We state in particular a Jordan decomposition of the evolution 
operator which is the key point in order to show that the fluctuations of the disordered 
Kuramoto model are not self-averaging. 



1. Introduction 

1.1. Synchronization of heterogeneous oscillators. Collective behavior of oscillators 
and synchronization phenomenon are the subject of a vast literature either in biological (neu- 
ronal models, collective behavior of insects, cells, etc.) or in physical contexts (see [25 11241 HQ] 
and references therein). While a precise description of each of the different instances in 
which synchronization emerges demands specific, possibly very complex models, the cele- 
brated Kuramoto model [U [39] has emerged as capturing some of the fundamental aspects of 
synchronization. 

The disordered Kuramoto model concerns a family of heterogeneous oscillators (or rotators) 
on the circle S := R/27rZ in a noisy mean-field interaction (that is the interaction is perturbed 
by a thermal noise), in the presence of disorder. More precisely, we make the assumption that 
the rotators are similar but not necessarily identical. In that sense, the model captures the 
fact that each rotator has a tendency to obey to its own natural frequency which may differ 
from one rotator to another. Those frequencies are chosen at random and independently for 
each rotator according to a probability distribution fi( dcj) on R; hence, this supplementary 
source of randomness will be considered as a disorder. One of the main characteristics of 
the Kuramoto model is that it presents a phase transition, as the coupling strength between 
the rotators increases, from an incoherent state where the rotators are not synchronized to a 
synchronous one where the phases of the rotators concentrate around a common value. 



1.1.1. The continuous model. The disordered Kuramoto model [251 [T], in the limit of an 
infinite population of rotators, is given by the following system of coupled nonlinear partial 
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differential equations (McKean-Vlasov equations): for all uj in Supp(/i) 

dtqt{e,oj)^^d^0qt{(^,u)-de(qt{e,Lo){{J*qt)^{9) + Lo)), t > 0, ^ e S, (1.1) 
with periodic boundary conditions and initial condition given by 

VL^eSupp(M), qt{e,Lo)de^j{de) (1.2) 

for some probability law 7 on the circle S. Here, 

{J * Qt)fj.iO) ^ -K sm{ip)qt{9-if,uj)dipfi{du}), (1.3) 

stands for the convolution of J(-) := -i^sin(-) with qt, averaged with respect to the disorder. 

Note that we are looking for solutions {t,6,u}) i-> qf(^6,uj) to (jl.ip - ()1.2p such that qt{-,uj) 
is a probability density for all fixed t and u: ■) ^ for all t > with j^qt{9,uj) dO - 1 
for all oj e Supp(/i). In that sense, the physical interpretation of (|l.ip ~ (|1.2p is the following: 
^ 6 S E [0, 27r] is the phase of the rotators, 7( dO) is the initial distribution of the rotators on 
S, -R' > is a nonnegative constant representing the coupling strength between the rotators, 
fjL is the probability distribution of the frequencies, and qt{9,uj) is the density at time t > of 
rotators with phase 6 and frequency uj. 

We stress the fact that evolution (jl.ip exhibits the following properties: it has quadratic 
nonlinearity, it is nonlocal, and if the support of /i is infinite, (jl.ip is in fact equivalent 
to an infinite number of coupled equations; indeed, in order to know q{-,uj) for one oj e 
Supp(//), it is necessary to know q{-,oj') for every uj' e Supp(/i). To fix ideas, one may 
for example consider the simple binary case where ^ - ^{5-i + 5i). Then (jl.ip reduces to 
two equations (one for +1, the other for -1) which are coupled via the averaged probability 
measure |(gt(^,+l) + qt{9 , -\)) dO . But for more general situations (/i = A/'(0,1) say) this 
would consist of an infinite number of coupled equations. 

Remark 1.1. It is immediate to see that (jl.ip exhibits the following symmetries: 

• Rotation invariance: ifqt{9,u}) solves (jl.ip so does qt{- + 6Q,uj) for any constant 9q e S, 

• Even symmetry: if qt{-9 , -uj)\t=o - Qti9,uj)\t=o, then it is true for all t> 0. 

In particular, the stationary solutions of (jl.ip will share these symmetries, see (jl.lSp . 

The Kuramoto model (jl.ip . as well as many generalizations, has been widely discussed 
mostly in the physical literature. We refer to [H [39] and references therein for a global 
review on the subject. We highlight in particular the active rotators model (see e.g. [37j ) 
which concerns similar interacting rotators in the presence of an additional potential. The 
dynamical properties for this model tend to be fairly more complicated than the model studied 
here. We refer to a recent result [19] that proves in particular the existence and stability of 
invariant normally hyperbolic manifolds for the active rotators model. 

Note that uniqueness of a solution to (jl.ip follows from standard arguments concerning 
fundamental solutions of parabolic equations ([21 [15]) and has been rigorously established 
in |17j § A]. Another proof of uniqueness can be found in [10] on the basis of heat kernel 
estimates under regularity assumptions on the initial condition. 

Existence of a solution to (jl.lj) may be seen as a consequence of the following probabilistic 
interpretation of (jl.ip . 
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1.1.2. The microscopic model. The system (jl.ip has a natural microscopic interpretation 
which is the fohowing: for all ^ 1, pick up independently at random a sequence of fre- 
quencies ui, . . . , LOjy e R according to the distribution /i and consider the following system 
of stochastic differential equations (whose solutions (%)i<jsg n are identified with their 
projection on S {9j mod 2tt)-^ < j < at)- 

K ^ 

yj^l,...,N, d9j,t = dt - - ^ sin(0j- 1 - 9i^t) dt + dBj^t , (1.4) 

i=i 

where at time t = 0, the rotators Oj are i.i.d. samples of the law 7 and {-Bj}j=i,...,Ar is a se- 
quence of N standard independent Brownian motions (which stands for the thermal noise in 
the system). We stress the fact that we have here two sources of randomness: the random- 
ness that chooses the frequencies (ujj)j^i and the randomness coming from the Brownian 
motions. Hence, one can see (II. 4p as a particular instance of interacting diffusions in a random 
environment. 

Evolution (jl.ip appears naturally as the large A^-limit of the system (II. 4p in the following 
way: let us define the empirical probability measure i^n of both rotators and frequencies: 

1 ^ 

t-i^TV,*:— E^(e,„a;,)eC([0,+oo),A^i(SxR)), (1.5) 

where S^g^^-j is the Dirac measure in {6,uj) and A^i(S x R) the set of probability measures on 
S X R. It can be shown (see [IHIISS]) that, under mild hypotheses, the sequence i^at converges 
in law (as a process), as goes to +00, to the deterministic limit t^vt such that 

I iyo{de,duj) = 7(d0)®/x(da;) 
where qt is solution of the McKean-Vlasov equation (jl.ip . 

Remark 1.2. We stress the fact that the above convergence holds for almost every choice of 
the initial sequence of frequencies > 1 (see j261 Th. 2. 5] J; due to the mean-field character 
of the diffusion (jl.4p . there is a self- averaging phenomenon.- the limit as N ^ +00 given by 
(jl.ip is deterministic and does not depend on the initial choice of the disorder {ujj)j ^ 1. 

The above convergence is a disordered generalization of known results about mean-field 
interacting diffusions. One can cite e.g. [16 j l21 j [22 l l29 j l32j for similar problematics concerning 
models without disorder. 



1.2. The fluctuation SPDE. In this paper, we investigate the asymptotic behavior as 
t +00 of the following stochastic partial differential equation: 

r/i = % + Lq.rjs ds + Wt, (1.7) 

where Lg^ is the linear operator arising from the linearization around the solution t >^ qt of 
nonlinear evolution (jl.ip : 

LgMO,io) ^d^eKO,io) - de {h{e,uj) ((J * qt)^{9)+uj) + qt{e,io){J * h)^{e)) , (1.8) 
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and where W is a Gaussian process, indexed by functions ip{9,uj) such that de(p{-,uj) e L^(S) 
for all u) 6 Supp(/x), with covariance 

Jf SAt r 
/ devi{e,u:)dev2{e,io)qu{e,u)deiM{du)du. (1.9) 

As we will see later in the introduction, the SPDE (jl.7p is the natural limit object in the 
Central Limit Theorem associated to the convergence of the empirical measure vn (II. 5p 
towards its limit v (II. 6p . as +00. Namely, the object of a previous work [26J was to prove 
that the fluctuation process 

^//V(^yiv-l/) (1.10) 

converges, in a weak sense, in an appropriate space of distributions on SxR, to the solution r/ 
of (ll.7p . The main motivation for this work was to study the influence of a typical choice of the 
disorder (ujj)j^i (quenched model) on the finite-size behavior of the microscopic Kuramoto 
model (jl.4p . As a matter of fact, it turns out that the limit process r] illustrates the dependence 
of the system (jl.4j) on the fluctuations of the frequencies ^i- as we will see, fluctuations 
of the disorder in (II. 4p compete with the fluctuations of the thermal noise and make the whole 
system rotate. What makes evolution (jl.7p relevant here is that its solution r] still captures 
this disorder-dependent rotation; at least numerically, one observes trajectories of the process 
r] that are compatible with the ones observed for the finite-size system ()1.4p (see Figures [2] 
and E]) . 

Hence, the main result of the paper (Theorem 12. 9p concerns the issue of the asymptotics of 
the solution ij of (jl.7p as time t goes to infinity. More precisely, the purpose of Theorem 12.91 
is to prove the following convergence, in an appropriate space of distribution, as t ^ +00: 

71+ in law 

f >V, (1.11) 

where the law of the limit v (that depends randomly on the disorder) is nontrivial and will 
be made explicit. In particular, this convergence might be an illustration of the fact that the 
self- averaging phenomenon for the law of large numbers observed in § 11.1.21 no longer holds 
at the scale of the central limit theorem: fluctuations in ()1.4p depend on the fluctuations 
of the disorder. Nevertheless, a rigorous derivation of the dependence in the disorder of the 
microscopic system (jl.4p would require to study this system for larger time scales (that is for 
time scales of order \/N instead of time scales of order 1, as it is done in this work). More 
precisely, one would need to study the joint limit of the empirical measure vn, as t and N 
goes to infinity. This has not been carried out her^ and would be a natural perspective for 
this work. 

The result we present here relies on a detailed spectral analysis of the unbounded evolution 
operator Lq defined in (II. Sp . using arguments from perturbation theory of self-adjoint opera- 
tors and of analytic semigroups. The main ingredient for this result is a precise understanding 
of the structure of the eigenspace of Lg associated to the eigenvalue 0, proving in particular 
the existence of a nontrivial Jordan block for the eigenvalue 0, relying on a priori estimates 
on the Dirichlet form associated to Lq and an extension of Lax-Milgram Theorem. 

Before going into the details of the asymptotics of (jl.7p . let us motivate our result in view 
of the recent advances made on the Kuramoto model. 



In this work, we first make N ^ 00 and then make f -> 00. 
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1.3. Motivations and problems. 

1.3.1. Notion of synchronization in the Kuramoto model. It is now well understood (see [36] ) 
that crucial features of evolution (jl.ip are captured by the order parameters ^ and ipt ^ S 
defined by: 

rte'^^^f e'^qt(9,Lo)d0dfi(uj). (1.12) 

JSxR 

The quantity rt captures in fact the degree of synchronization of a solution (the profile % = ^ 
for example corresponds to rt = and represents a total lack of synchronization) and 'ipt 
identifies the center of synchronization: this is true and rather intuitive for unimodal profiles. 
In this framework, it turns out that synchronization reads in the existence of nontrivial 
stationary solutions q to (|l.ip . which can be parameterized by the order parameter r defined 
in (jl.l2p : namely, following [36] (see also [111 p. 117]), if /i is symmetric, any equilibrium in 
(11. ip may be written as q{6 + 9q,uj) for any fixed arbitrary 0o ^ S where 

S{e,u:,2Kr) 
'^^^'"^^ Z{u:,2Kr) ' 

for 

5(0,w,x) :=e^(^'^'^')[(l-e^'^^) T^-^^"''^-^) dn + e^^"^ T^^e-^^^^dJ, (1.14) 
V Jo J 

where G{u, u},x) - x cos(u) + 2tjn, Z{oj, x) - /g S{9, uj, x) dO a normalization constant in order 
to insure that q is indeed a probability density on S, and r e [0, 1] satisfies (jl.l2p with q given 
by ^M- 

, , r /■^cos(0)5(e,w,x)d0 , ^ , ^ 

r^^^{2Kr), where ^'^.(x) / ^-^ \!, \ /"(do;)- (1-15) 

Discriminating between the possible solutions r e [0, 1] to the fixed-point relation (|1.15p . one 
has two possibilities: 

• The case r = is always a solution to (jl.lSp and corresponds to the constant density 

• Any solution q with r > is called a synchronized (nontrivial) solution. Any easy 
calculation of the derivative of ^'^(•) at shows that such solutions exist at least 

when the coupling strength K is greater than K i+tui^ ) ' ti^s^t case, due to 

the rotation invariance (Remark ll.ip . each solution r > of ()1.15p corresponds to a 
whole circle of synchronized stationary solutions 

Msync---{qi- + Oo,u;)-eo€S}. (1.16) 

1.3.2. The case with no disorder. In order to motivate our results, let us give a brief review 
of the case where there is no disorder, that is /i = (Jq- In that case, (jl.ip reduces to: 

dtqt{9) - Id'^MG) - de iqt{e){J * %)(^)) • (1-17) 

As well as in the disordered case, any non-trivial stationary profile can be written as go(^ + ^o) 
where 

g2_ft'ro cos(e) g2A'rocos(e) 
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where {K > 1) ro is the only nonnegative solution of 



ro = ^'o(2i^ro), where ^'0(2;) := 



/g cos(^)e' 



(1.19) 



In that case, since ^'o is strictly concave ([34, Lem. 4]) and 9,.g^'o(2i^ro) = K, the phase 
transition is obvious: for ^ 1, tq = is the only solution to (jl.lOp and ^ is the only 
stationary solution whereas for K > 1 this solution coexists with a unique (up to rotation) 
synchronized solution (corresponding to the unique > solution to (11.191) ). 

Recent advances have been made concerning the dynamical properties of the non-disordered 
Kuramoto model. In particular in [5], the linear stability of the synchronized solution (jl.lSp 
is proved, with explicit spectral gap estimates. We will rely on these results for the rest of 
this paper and we refer to § l2.3l for more details. What is more, the whole long-time dynamics 
of (ll.lTh is addressed in [18] , discussing in particular the existence and the nature of a global 
attractor in both cases K ^ I and K > 1 and proving non-linear stability of synchronized 
solutions for K > 1. 

1.3.3. Perturbing the non disordered case. A natural issue is then to understand in what way 
adding the disorder (wj)-,- ^ 1 may perturb the results mentioned for the non-disordered case. 
In this direction, the question of linear and nonlinear stability of synchronized solutions in 
the disordered case is addressed in p/7j, at least in the case where the disorder is small. In 
that extent, one can distinguish two regimes: 

• The case fi not symmetric: from an intuitive point of view, in this case, there are more 
rotators that are likely to rotate in one given direction than in the other direction; 
the whole system has a tendency to rotate in this direction. This phenomenon can 
be made rigorous and reads in fact on the continuous model (jl.ip : for small disorder, 
the existence of stable rotating- waves which solve (jl.ip is proved in [T7j, with explicit 
computation of the speed of rotation. 

• The case fi symmetric: in this case, the speed of rotation of the rotating waves in 
the continuous model (jl.ip is equally zero. Nevertheless, one can prove (^17j) by 
perturbations arguments the linear and nonlinear stability of synchronized stationary 
solutions, once again for small disorder (with explicit estimates on the size of the 
disorder). 

1.3.4. Finite system: disorder-induced rotation and non self-averaging phenomenon. As al- 
ready mentioned, the study of the SPDE ()1.7p is initially motivated by the Central Limit 
Theorem associated to the convergence of the empirical measure defined in (|1.5p and 
is linked with the behavior of the microscopic system ()1.4p for large but finite N. Indeed, 
Balmforth and Sassi [3] pointed out on the basis of numerical computations the following 
remarkable phenomenon: even if the distribution fi is symmetric, the fiuctuations of a fixed 
chosen sample of the disorder (wi,... ,ujn) makes it not symmetric and induces a slow ro- 
tation of the center of the synchronization in one direction at a speed of order N'^l"^ that 
depends randomly on the sample of the disorder (Fig. [T]and[2]). To fix ideas, one may think 
of the simple case where the distribution of the disorder is binary: /i = ^ (<^-i + In a 
finite sample (wi, . . . , wat) of spins in {±1}^, there may be a majority of positive spins (equal 
to +1) with respect to the negative spins (equal to -1), so that the rotators with positive 
frequency have a tendency to induce a global rotation in their sense. This non self-averaging 
phenomenon at the scale of fiuctuations (to compare with Remark II. 2p can be noticed by 
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Figure 1. We plot here the evolution of the marginal on S of z/^r for = 600 
oscillators {fi - ^{S-i +Si), K - 6). The oscillators are initially chosen inde- 
pendently and uniformly on [0, 27r] independently of the disorder. First the 
dynamics leads to synchronization of the oscillators (t = 6) to a profile which 
is close to a nontrivial stationary solution of McKean-Vlasov equation (jl.ip . 
Secondly, we observe that the center ipN,t of this density moves to the right 
with an approximately constant speed; what is more, this speed of fluctuation 
turns out to be sample-dependent (see Fig. [2]). 



computing the finite-size order parameters {rN,tT'4'N,t) (that is the microscopic equivalents of 
the macroscopic order parameters {rt,ij)t) defined in (jl.l2p ): 

1 ^ r 

where is the empirical measure of the system defined in (jl.Sp . In particular, the order 
parameter ^l^N,t captures the position of the center of synchronization for finite (see Fig- 
ure [1]). It is seen on Figure [2] that 1 1-^ V'A'.t has an approximately linear behavior whose slope 
depends on the choice of the disorder. 

1.3.5. Central Limit Theorem. The first step to tackle this phenomenon was to prove (see [26j) 
that the following convergence holds 

VN V (1.21) 

in the space of continuous processes with values in iS'(S x R) where rj is the solution to (jl.7p 

r/i = r?o + ^ Lg^7]s ds + Wt, 

and where 

• 5'(S X R) is understood as the usual space of tempered distributions on S x R, 

• W is the Gaussian process whose covariance is given in (|1.9p . 

• Lq^ is the operator defined in p.Sp . 
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Figure 2. Trajectories of the center of synchronization for different re- 
ahzations of the disorder (// = |((5-o.5 + ^0.5), K = 4, N - 400). We observe 
here the non self-averaging phenomenon: direction and speed of the center 
depend on the choice of the initial A'^-sample of the disorder. Moreover, these 
simulations are compatible with speeds of order N^^^"^. 



• the initial condition r/o is independent with W. 

It is important to note that the convergence that is considered in (jl.2ip is more elaborate 
than a usual convergence in law, since this convergence is understood in a quenched set-up, 
that is for fixed choice of frequencies > 1. In fact, the precise notion of convergence used 
in [26] for (jl.2ip is not really relevant for the purpose of this paper. The interested reader 
may find more details in [2_51 Th. 2.10]. What we only need to retain here is that, in view of 
the fact that there are two sources of randomness in the microscopic model (jl.4p , the process 
T] also depends on two sources of randomness at the limit. More precisely, the limit object 
rj captures the dependence of the disorder through its nontrivial mean-value: there exists a 
Gaussian process u 1-^ C{uj) with covariance 

V(/Ji,(/72 : S X R ^ R, rc{ipi,(p2) ^ CoY^^J^ipi{-,u;)d'y, J^ip2{-,uj)djj (1.22) 

such that for fixed u, the initial condition rjQ = rj^ in (jl.7p may be written as 

r,^^X + C{uj), (1.23) 

where X is a centered Gaussian process with explicit covariance, and C(uj) is the constant 
mean- value of r]Q . Let us give a microscopic intuition of the decomposition p.23p : at time 



LARGE TIME ASYMPTOTICS FOR THE FLUCTUATION SPDE IN THE KURAMOTO MODEL 



9 



t - 0, for N ^ 1 and (/? : S x R -j- R, rjiy Q^cp) defined in (jl.lOp may be decomposed into 
W,o('/') = i E ^i) - £ „ 9'(^,^)7(d6')^(da;)], 

= E (^(% , ) - ^ 9^(0, )7( d^) ) 

+ £ ^,)7( d^) - (^(^,^)7( de)/^( d^)), (1.24) 

In particular, by definition of X^, this process is centered for all choice of (ujj)j^i and 
converges to some X for almost every {ujj)j^i ([6l Th. 27.3]). As far as the process Cn 
is concerned, it is easily seen that it converges in law (w.r.t. the disorder) to the process C 
with covariance given by ()1.22p . What is important to understand here is that the random 
variable C(-) precisely captures the fluctuations of the disorder. For example, let us consider 
the case where ^ = ^ i^-uo + <5wo)' some wq > 0. As we will see in Theorem 12.91 the relevant 
quantity is in fact the process C+ (indexed by functions "0 • S R) that is the restriction of 
the process C to the component on +ijJo: 

VV:S-R, C+,^:=C^l_^. (1.25) 

Thanks to (11.240 . C+ is the limit in law of the microscopic process Cm.+ (indexed by functions 
: S R) defined by 

V*. C„,W.(X*(.)do)^|(l,_,,-l):=(XW)d7)^. (1.26) 

Here, a at is exactly the (centered) number of frequencies among (oji, . . . jLUn) that are positive, 
so that Cn^^ captures the lack of symmetry of the initial chosen disorder: if oat = 0, there are 
as many positive frequencies as negative frequencie^, whereas the case on > (resp. <0) 
represents the case of an asymmetry in favor of positive (resp. negative) frequencies. 

Numerical simulations (Figure [3|) of the process given by (|1.7p show a non self-averaging 
phenomenon for the fluctuation process compatible with the one observed for the rotation 
of the finite system (Figure [2|): as seen in Figure [3l the trajectories of r/(sin) have an ap- 
proximately linear behavior, which slope depend on the realization of the initial mean-value 
C. Hence, a way to understand the dependence of the microscopic system (11.40 in the initial 
asymmetry of the disorder is to analyze the dependence of the fluctuation process ij (11.70 in 
the mean- value C of its initial condition. More explicitly, the key point of this paper is to un- 
derstand how different initial conditions in evolution (11.70 may lead to distinct approximately 
linear trajectories of the fluctuation process. 

We tackle this issue via a spectral analysis of the evolution operator Lg^ defined in (II. Sp via 
arguments from the general theory of perturbations of self-adjoint operators ([23]) and analytic 
semigroup of operators ([271 [33]) aiid usual techniques about SPDEs in Hilbert spaces ([9j), 
using the precise knowledge we have about the same linear operator Lg^ in the case without 
disorder (see [5]). 



'in such case, the number A'' of rotators is necessarily even. 
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77t(sin) 




Figure 3. We plot here the evolution of the process rjt{sm), for different 
realizations of C; the trajectories are sample-dependent and compatible with 
the ones observed in Fig [2j 



1.4. Organization of the paper. The paper is organized as follows: in Section [2l we precise 
the main set-up for the study of the SPDE ()1.7p and state the main results. In particular, 
Theorem 12.51 and Theorem 12.71 deal with the spectral properties of the evolution operator 
at least when the disorder is small. Secondly, we state the main result of this paper: 
Theorem 12.91 establishes the linear asymptotics of the fluctuation process solution of ()1.7p . 

Section [3] is devoted to prove Theorem 12.51 In Section [5] we prove Theorem 12.71 whereas 
the main result of the paper. Theorem 12.91 is proved in Section [5l 

2. Main definitions and results 

In what follows, the following standard notations will be used: for an operator F, we 
will denote by p{F) the set of all complex numbers A for which A - -F is invertible, and by 
R{X,F) ■■- (A-F)~ , A e piF) the resolvent of F. The spectrum of F will be denoted as 

2.1. The evolution operator Lq. The main issue of this paper is to prove asymptotic 
results for the solution of the SPDE (jl.7p as t ^ +oo; this fact is deeply linked to the spectral 
properties of the operator Lg defined in (jl.Sp . We will restrict ourselves to the stationary 
case, that is when q\t=o - It is equal to the synchronized (non-trivial) stationary solution q 
(jl.lSp of evolution (jl.ip . In this case, the object of interest is the stationary version of p.Sp . 
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that is 

Lh{e,io) ^d^eKG:^) - de {h{9,uj) ((J * q)^{0) + io) + q{9,io){J * h)^{e)) . (2.1) 
The domain T> of the operator L is given by: 

|/i(6l,a;); Va;,6l H> /i(6»,tj) eC^(S), [ /i(6l,a;) d6l^( dw) = ol . (2.2) 

I. JSxR J 

Remark 2.1. We point out that the choice of the domain T) of L is of prime importance 
for our study of evolution (jl.7p . In another context, one encounters the same operator L in 
the study of the linear stability of the stationary solution q since the linearized evolution of 
(jl.ip around q is precisely given by dthf - Lhf. The natural domain for this latter evolution 
(see [T7]j is then 

h{0,uj);yuj,e^h{e,uj) eC^{S), Vw, J^h{e,uj) dO^Oy (2.3) 

Indeed for all lo, q{-,u}) is a probability density on S so that perturbing by elements of domain 
()2.3p enables to remain within the set of functions with integral 1 on S. 

But here, the key point is to understand that evolution (jl.7p does not live in domain ()2.3p 
since the fluctuations of the disorder are precisely captured by the fact that r] has a non-trivial 
mean-value C(uj) for fixed lo. We will henceforth work with the domain T> defined in ()2.2p 
instead of ()2.3p . This choice of domain (apart from being technically more demanding) is all 
the more relevant since we will see that the non self- averaging phenomenon holds in ()2.2p and 
not in (j2.3p (see Remark \2.6\ below). 

For the rest of this paper, we fix > 1 and we restrict ourselves to the case where 



where wq > is a fixed parameter. This assumption (j2.4p on fi appears to be quite restrictive, 
but generalizing parts of the results we present here to more general distributions fj, does not 
seem to be straightforward. We refer to § 12.71 for a discussion on this topic. 

The first main goal of this paper is to state a spectral decomposition of the operator L 
defined in (12. ID . based on perturbation arguments from the non-disordered case fj, - 6q (see 
5 10:2] and §[23]). 



2.2. Distribution spaces. The spectral analysis of the operator L will be mostly carried out 
in spaces of distribution that have H^^ regularity w.r.t. 6. But the precise study of L requires 
to introduce weighted version of H'^ that we define here. We first focus on weighted-Sobolev 
spaces of functions 9 ^ h{9) on S (§ 12.2. ip and then introduce the corresponding spaces for 
functions with disorder {6,Ijj) ^ h{6,io) on S x Supp(;u) (§ I2.2.2P : 

2.2.1. Weighted Sobolev spaces: For any bounded positive weight function fc(-) on S such that 
/g k{9) d9 = 1, we may consider the space closure of C(S) w.r.t. the norm: 

The decomposition of h into the sum of Span(A;) and its orthogonal supplementary in may 
be written as: 

h 



= (_^/i)-A: + /io, (2.6) 
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where /g ho - 0. Since /iq is with zero mean value, each of its primitives are 27r-periodic. In 
particular, we can consider H^^ the closure of C(S) with respect to the following weighted 
Sobolev norm: 

1 

m-l,k:-(l fh]\ L^]\ (2-7) 



where Hq is the primitive of /iq on S such that /s = 0- Note that one can understand the 
spaces H^^ as part of a Gelfand-triple construction (see Appendix |A] for a precise definition). 
In particular, we will make a constant use of the space (that is for k{-) - qo{-) where qo 
is the stationary solution (jl.lSp of the non-disordered system) which is the natural space (see 
Prop. 12. 4p for the study on the Kuramoto operator Lq^ (|2.11[) in the non-disordered case. 

Remark 2.2. In the case of a constant weight k{-) = we will write (L^, || • II2) and (H""*^, || • 
II _i) instead 0/ (L^, II "11 2 j_) and (H^, || ■ ||_-|^ j_). 

27r '2ir 2n ' ^Ti 

2.2.2. Weighted Sobolev spaces (with disorder): The natural space in where to study the 
operator L is the space of functions h in V such that each component /i(-, w) live in a certain 
^k{- uj) ^ weight k{-,io) (which may depends on a; e Supp(/x)). More precisely, for any 
family of positive weight functions , ui)) ^^g^^^^^^-^ , we denote as the closure of P w.r.t. 
the norm: 



II fJ.,-l,k 



We will also consider the analogous averaged weighted L^-spaces, that is, for any family of 
positive weights (A;(-,a;))(^£Supp(^)) the space ^ given by the norm: 

Remark 2.3. In the particular case of k{-,uj) = ^ for all uj e Supp(^), we will write 

H,/ := H,/ 1 



2 



and the corresponding norm will he denoted as || • || jj . We will also write (L^, || • ||^ 2) instead 

The main theorem concerning the operator L (Theorem 12. 7p will be stated in for the 
ease of exposition but its proof will require the introduction of weighted-Sobolev spaces H~^, 
for nontrivial weights k. 

2.3. The non-disordered case. In the context of the Kuramoto model without disorder, 
the linearized operator Lq^ around stationary solution qq (see § I1.3.2P , with domain 



(S);^n = o} (2.10) 



is: 



1 



LqoU := -dgu - de [qo{J * u) + u{J * qo)] . (2.11) 



2 

In [3], it is mainly proved that Lq^ is essentially self-adjoint in H~^: 
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Proposition 2.4 (Bertini-Giacomin-Pakdaman). (Lqg,Po) ^-^ essentially self-adjoint inJl^^. 
The spectrum of (the self-adjoint extension of) Lg^ is pure point lying in (-00, 0); is in 
the spectrum, with one- dimensional eigenspace (spanned by dgq^). Moreover, the distance 
XxiLgg) between the eigenvalue and the rest of the spectrum is strictly positive. 

Proposition 12.41 is crucial for our study of operator L (j2.ip . We will refer to it when required. 

2.4. Non self-averaging phenomenon for the operator L and existence of a Jordan 

block. The linear trajectories that depend on the initial condition as observed in Figure [3] 
are reminiscent of an analogous deterministic finite-dimensional example: let us consider the 
2-dimensional evolution - -^(y(t))' ~ (oo)- Then it is trivial to see that the 

solutions of this system are linear in time: yo as t ^ 00. This can be easily generalized 

for matrices L n^n with spectrum in {A e C; 9^(A) ^ 0} such that the Jordan decomposition 
of L admits a ( [] q )-block; this fact is precisely equivalent to the existence of x and y such 
that Lx - and Ly - x. The purpose of the first main theorem of this paper is to prove an 
analogous existence of a Jordan block for the operator L: 

Theorem 2.5. For any fixed ujq > 0, if q is the stationary solution in (ll.lSp . then 

Ldeq^O. (2.12) 

Moreover, there exists p^D such that 

6 S, Vtj e Supp(/i), Lp{e,u) ^ deqie,u;). (2.13) 
In particular, the characteristic space of L in is at least of dimension 2. 

Remark 2.6. Equality (|2.12p is a direct consequence of the rotation invariance in (|1.4p 
( Remark \ 1.1]} . Note also that p{-,u}) found in (|2.13p is with non-trivial mean value for all 
oj e Supp(^). Indeed, we believe that J^p{-,uj) is in fact equal to this fact is derived from 
non-rigorous computations and verified by numerical simulations. In other terms, if we had 
chosen domain ()2.3p instead of domain (j2.2p . we would not have found such p; the Jordan 
block (00) the matrix representation ()2.17p of the operator L does not exist on the domain 



Theorem 12.51 is proved in Section [3l 

2.5. Spectral properties of L and position of the spectrum. The second goal of this 
paper is to prove that L generates an analytic semi-group of operators with spectrum confined 
in the part of the complex plane with negative real part: 

Theorem 2.7. In the Hilbert space H^"*^ defined in Remark ] 2. 31 the operator {L,T>) is densely 
defined, closable, its closed extension having compact resolvent. In particular, its spectrum 
consists of isolated eigenvalues with finite multiplicities. 

Moreover, for all K > 1, for all a e (0, |), for all p e (0, 1), there exists cj* = uJi,{K, a,p) > 
such that, for all < ojq < UJ^,, the following is true: 

• The spectrum of L lies in a cone Ca with vertex and angle a 

C« |a 6 C; I + a ^ arg(A) ^ y - a| c {z e C; lH(z) ^ 0} ; (2.14) 

• There exists a' e (0, ^) such that L is the infinitesimal generator of an analytic semi- 
group defined on a sector Aa' ■- {A e C, | arg(A)| < a'}; 
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• the dimension of the characteristic space in is exactly 2, spanned by dgq and p, 
where p is defined in Theorem \2.5\ 

• the eigenvalue is separated from the rest of the spectrum at a distance Xk{L) - 
X{L,K,p) at least equal to p • min (Ai^(LgQ), ^e"^^''°), where Lq^ and ro are defined 
in §{MM 

Note that Theorem 12.71 relies on perturbation arguments of the non-disordered case men- 
tioned in S I2.3t in that sense, the spectral gap Xk{L) found in Theorem 12.71 depends on the 
spectral gap Ai<:(Lg„) for the non-disordered case. 

As a consequence of Theorem 12. 7| there exists a decomposition of the space into the 
direct sum 



(2.15) 



where Go is of dimension 2 (spanned by dgq and p) such that the restriction of the opera- 
tor L to Go has spectrum {0} and the restriction of L to G<o has spectrum cr(L) \ {0} £ 
{A e C; 9^(A) < 0}. We will denote as Po the corresponding projection on Go along to G<o , 
and E:o - I - Pq- In particular, there exist unique continuous linear forms Ig^g and £p such 
that for all h e H;;^ 

Poh:^ldeq{h)deq + lp{h)p. (2.16) 
To fix ideas, one may think of the following infinite matrix representation for the operator L: 

I 



L = 











PqLP. 



<0 



Pia LPi 



<o 



\ 


/ 


1 


edoALP<^) ] 


}deq 


|go 










••• 


}p 


























' G<o 


/ 


\ 





1 







(2.17) 



Span(<9eg,p) 



G<0 



Note that the second line in the matrix representation (j2.17p of L is indeed equally zero since 
for all h e H~^, Lh is of zero mean value on S; in particular ip{Lh) - for all h e H^^. 

Remark 2.8. Any element h - (^(^) '^))ees u;eSupp(/x) '-^^ identified in our binary case 
(|2.4p with a couple (/i+(^), /i-(6'))0es- Moreover, any h e H^^ can be decomposed according to 

h = idgq{h)dgq + £p{h)p + P^q h. 

Let us integrate the latter decomposition w.r.t. 9. Since j^Lu - for all u € D, we have 
J^P-qH- so that one can actually find an explicit formulation for the functional ip: 

W -_k!!i-_k!L. (2.18) 

JsP+ JsP- 

The last equality in ()2.18p is due to the fact that /g (h+ + /i_) = /g {p+ + p-) - 0. 
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2.6. Long time evolution of the fluctuation SPDE. We now turn to the main result 
of the paper, which concerns the asymptotic behavior of the fluctuation process r] defined in 



where p+{0) p(0,ujo) is defined by (|2.13|) . Moreover, in view of the structure of the initial 
condition r/o (recall (|1.23p ). the speed v(-) is a Gaussian random variable in uj with variance 



Hence, Theorem 12.91 gives an explicit expression of the speed of rotation of the fluctuation 
process. 

2.7. Perspectives. 

2.7.1. Possible generalizations. In [4, § 10.2, p. 47], Balmforth and Sassi found on the basis 
of numerical experiments that if we get rid artificially of the asymmetry between positive 
frequencies and negative frequencies, the disorder-induced rotation in (jl.4p disappears. More 
precisely, the simulations carried out in [3] are ruled by the following principle: consider a 
sample of 2N rotators (A^ ^ 1) and pick up at random the frequencies for exactly 
rotators. For each of those N frequencies Wj, attribute to one the remaining N rotators the 
exact opposite frequency -Wj. The conclusion of ^ is that for this choice of frequencies, there 
is no rotation, no matter how the first frequencies are sampled. We stress the fact that 
we retrieve this phenomenon in Theorem 12.91 in the case where the first N rotators are i.i.d. 
samples of /i = | (5-a;o + ^uio)^ since in that case the quantity a2N defined in ()1.26p is equally 
zero for all ^ 1 and the consequent limit speed v is equally zero in this case. 

One could hope to generalize the results of the paper in at least two directions. Firstly, we 
restrict ourselves to the binary case /x = \{5-uq + ^ujq)- Note that the proof of Theorem 12.51 
concerning the existence of a Jordan block (although written in this particular case for the 
reader's convenience) is not specific to this restrictive case: one could easily rewrite the same 
proof in the case of more general distributions /x (even with unbounded support), satisfying 
appropriate integrability conditions in and in oo. 

The main restriction on ^ concerns Theorem 12.71 the hypothesis /i = +'^wo) is critical 

for its proof. In that sense, there are some deeper considerations which may prove possible 
generalizations of Theorem 12.71 to be difficult: indeed, the key argument of the proof is based 
on the fact that perturbing a finite dimensional kernel of an operator ^ by a sufficiently small 
perturbation B leads to a kernel for the operator A + B with the same finite dimension. But 
for distributions more general than binary distributions |(5-(^o +5wo)' kernel of L is likely 
to become of infinite dimension, so that those perturbation arguments cannot be applied. 

Secondly, Theorem 12.71 is only proved for small disorder whereas one would expect it 
to be true even for large disorder. It is indeed natural to believe that the non self-averaging 
phenomenon seen in Figure [2] not only holds for large disorder but would even be more 
noticeable in that case. However, since Theorem 12.71 relies on perturbation arguments, to 
prove similar results for large ujq would seem to require alternative methods. 




^^^^^deq-^v{uj)deq, ast^+oo, (2.19) 





(2.20) 



16 



ERIC LUgONi 



2.7.2. Other possible time scales. More importantly, Theorem 12.91 gives an expression of the 
speed at which the hmit of the fluctuation process ij rotates. Since in this work, we are only 
considering time scales of order 0(1) (that is for arbitrary but fixed time horizon T), this 
may not be the right speed of rotation for the finite-size system: on larger time-scales (or 
order 0{\/N) say), there may be accumulation of fluctuations which may lead to a different 
expression of the speed of rotation. Hence, to understand this phenomenon properly, an 
interesting result would be the study of the rescaled empirical process {u ^^)t,N, as and 
t are large. This has not been carried out in this article and would be a natural perspective 
for this work. 

3. On the existence of a Jordan block for L (Proof of Theorem 12.51) 

The purpose of this section is to prove Theorem 12. 5[ that is the fact that the operator L 
defined in (|2.ip has a Jordan block of size at least 2. The symmetry of the system (recall 
Remark I l.ip leads to consider the set of distributions which are odd w.r.t. the variable {0,ljj) e 
S X Supp(/x): 

O := {h; V(^,w) e S x Supp(^), h{-9,-uj) = -h{9,uj)} . (3.1) 
We also denote by M the set of functions that with zero mean-value for all uj e Supp(//) (recall 
the definition of V in STlh): 



A/':=|/ieP; Va;eSupp(^),^/i(6l,w)d6l = o|. (3.2) 

In the following straightforward lemma, whose proof is left to the reader, we sum-up the basic 
properties of the stationary solution q (recall (jl.l3p ): 

Lemma 3.1. The following statements are true: 

(1) dgq^O, 

(2) Ifh^O then Lh e O, 

(3) dgq e M, and for all h^V, LheJ\f, 

(4) There exist some positive constants c and C such that for all 6 € S, lo € Supp(/i), 
0<ci^q{e,uj) ^C, 

(5) For all 9 € S, CO € Supp(/u), 

^deq{9, u) = q{9,u) ((J * q)^ + u) + k{u;) , (3.3) 

where - ^^§(jjy, and Z{uj) - Z(uj,2Kr) is the constant of normalization defined 

in (fLTSj) . 

Note in particular that the fact that d^q e O can be seen as a consequence of Remark II. li A 
direct calculation easily shows that dgq is in the kernel of L (it corresponds to the rotation 
invariance of the problem). The rest of this section is devoted to prove the existence of an 
element p &T> such that Lp - dgq. 

We recall here the definition of the weighted Sobolev spaces introduced in § I2.2.2t we make 
use here of the space (H^^^, (■ , ■)^ ^ ^) defined in (j2.8p in the case of k - q and of the space 

(L^, (• , •)^ 2) defined in Remark 12.31 The main result is the following: 

Proposition 3.2. For every luq > 0, in the binary case (j2.4p . for every v e H^^ n O (and in 
particular for v - dgq), there exists some p e n O such that 

yi 6 H^;^ n O, {Lp , = {v , 0^,-1,, • (3.4) 
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Moreover, in the case v - dgq, any p that satisfies ()3.4p is in fact a regular function (p(-,uj) e 
C°°(S) for all to e Supp(/i) ) and is a classical solution to (j2.13p . 

Remark 3.3. The scope of Proposition \3.S\ is more general than the restrictive case of a 
binary law = ^ ('^-i^o + ^i^o)> indeed the following proof works for more general distributions 
H, the only additional requirement being integrability condition^ in and +oo, see Remark \3.9l 

Proof of Proposition 13.21 relies on several lemmas: 

Lemma 3.4. For h € O nV, I € O nV, let us introduce the Dirichlet form 

£L{h,l):^{Lh,l)^^_,^g. (3.5) 
£l{; •) is well defined on V{£l) := (L^ n O) x (H^i^ n O) and one can decompose £"/,(•,•) into: 
V(/i, e V{£l), £L{h, I) = T{h, I) + Kl{h) ■ £{l), (3.6) 
where T{-, •), bilinear form on T>{£i,) and £{■) linear form on H^^'^ n O, are defined as follows: 

y{h,l)eV{£L), r{h,l):^-l f -dXdfi+ f K{-)^dXdfi, (3.7) 

2 JSxR q JsxR 

V/ e U-^]g nO, ^iO- sin(-) dA dfi, (3.8) 

where k and L in (13. 7p are respectively defined in (13. 3p and as the primitive of I e H^^^ such 
that /g = for all oj e Supp(^) (recall §\EEM)- 

Lemma 3.5. For all continuous linear form f on H^^'^ n O, there exists some pi e L^ n O 
such that for I e ^ n O 

r(pi,o = /(0- (3.9) 

Lemma 3.6. The linear form i(-) defined in (j3.8p can be expressed as a scalar product on 
Hjj^g n O: there exists p2^T>nO, for all I e H^^^ n O 

£{l)^{Lp2,l)^^_,^g. (3.10) 

Let us admit for a moment Lemmas 13.41 13.51 and 13.61 and let us prove Proposition 13. 2t 

Proof of Proposition \3.^ Let u be a fixed element of H^^^ n O. Applying Lemma 13.51 to the 
continuous linear form /(/) = {v , I) ^ _^ g, there exists some pi e L^ n O such that T{pi,l) - 
{v , l)^_iq, which gives using Lemma 13.41 and Lemma 13.61 

0^.,-i,g = r(pi,/), 

We can conclude that the variational formula (|3.4p is verified for the following choice of p: 

p:^Pi-Ki{pi)p2€LlnO. (3.11) 

Let us prove now that such p is in fact a regular function in 9: since p2 ^"D is regular in 0, it 
suffices to prove that for all uj e Supp(/i), 6 >->■ pi{6,uj) is (in fact C°°) in 6. We start from 
the definition of pi : 



those conditions are obviously satisfied in the binary case (|2.4[) . 
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Since this true for all I e H^^^ n O, thanks to the expression of F in ()3.7p . we obtain that for 
any fixed u e Supp(/i), for Lebesgue-almost every 6 e S: 



\Jo q(u,ijj)'^ I Jo q(u,ijj) 



4 q{0,uj) \Jo q{u,uj)^ J Jo q{u,u}) 

where Q{-,uj) is the primitive of dgq{-,uj) such that /g ^^r^ - 0- But, using that q is bounded 
and C°° in 9 and that pi(-,uj) e L^, the primitive /o^ gli^tl))^ ^ ^ version. So, thanks 
to (j3.12p . pi has a version. So, the right-hand side of (j3.12p is a least C^, and so does 
Pi. The same repeated argument shows that pi is C°° in 9. That concludes the proof of 
Proposition 13.21 □ 

It now remains to prove the three lemmas: 

Proof of Lemma \3m Let us prove equality (|3.6p : since £ is a primitive of I, one has 

1 r {deh)C _ 1 rhl_ 1 f h£r. 

Using p.3p . for a; e Supp(/u) 

1 r{deK)C I r hi fhC..^ , , t ^ f 

2 7s^^-2 7s7"7sT^^^^'^^-^-^""^"'^^"^7s^- 

Thanks to the expression of Lh in (|2.ip . we obtain 

£L{.h,l)^-\f -+ f k{-)\c- f {J*h),£, (3.13) 

2 JSxR q JSxR JSxR 

= r(/i,/)- f {J*h)^c. 

JSxR 

Lastly, integrating by parts the last term in (13.13P and expanding the cosine function (recall 
J(-) = -Ksin(-)), we obtain: 

^I jJ . k),C = K(/^^^c„B(.)MAd,) (^^^co,(.)/>dAd,) 

+ k( [ sin(-)/dAd^U [ sin(-)/idAd/iV 

VJSxR / VJSxR / 

But, since / e O, the first term in the latter expression is zero. The result (13. 6p follows. □ 



Proof of Lemma \3.5[ In this proof, we use the following extension to Lax-Milgram Theorem 
(cf. [3S chap. Ill]): 

Proposition 3.7 (Lions-Lax-Milgram) . Let {Ti, |.|} be a Hilbert space and {Q, || . ||} a normed 
linear space. Suppose F : x ^ ^ R is bilinear and that F(-, if) is continuous for each ip € Q. 
If there exists some constant C > such that 

inf sup |F(/i, (/?)| ^ C > 0, (weak coercivity), (3-14) 
llvlhil/il s: 1 

Then for each f ^ G' there exists some p ^T-L such that T{p,Lp) - f{(p) for all (p ^ Q. 



LARGE TIME ASYMPTOTICS FOR THE FLUCTUATION SPDE IN THE KURAMOTO MODEL 19 

The principle of the proof of Lemma 13.51 is to show that the bihnear function F defined in 
(13. 7p satisfies Proposition 13.71 for 

"Hi^L^nO, endowed with ||-||^2' (3.15) 
g:=H^]gnOnL°°(SxR), endowed with Ml^,-i,q- (3-16) 

Namely, we have the following: 

(1) For each / e Q, r(-,/)is continuous on n O: indeed, for the first term of T{h,l), we 
have 

f —dXdfi ^CIUL f l^|dAd//^C||/i||„2. 

JSxR q JSxR "^'"^ 

And for the second term, using the boundedness of q: 

f \k(-)-!^C dXdfi^C [ |/i/:|dAd^^C||/i||„2|U||„ i„. 

(2) r is weakly coercive: let us fix / e ^ such that IM ||^ _i ^ = 1- 
Let us choose h = gC € hj^ n O, where for all uj e Supp(/u), g{-,uj) is a 27r-periodic 

function to be defined later. Then, by integration by parts in the equality (j3.7p 

Consider now for fixed u: e Supp(/^t) the following first order ODE, with periodic 
boundary condition: 

W(-,w) + Ac(a;)4^ = ^^, with/(0,L^) = /(27r,w). (3.18) 

Then for any u e Supp(/i) \ {0}, an explicit calculation (left to the reader) shows that 
there exists a unique solution to (j3.18p . 6 ^ f(6,u}). 

Remark 3.8. In the case oj - 0, (j3.18p reduces to \def - which is incompatible with 
the condition /(O) = /(27r), since J^^dO > 0: there is no such 2iT-periodic solution 
in the case w = 0. 

Moreover, it is straightforward to see that || |/(-,cj)| d/x ||^ ^ ^ C, for some con- 
stant C > 0. 

Remark 3.9. It is easy to see that f{-,oj) is not bounded as uj ^ and oo +oo; 
thus, for general distributions fj,, the same control on f requires additional integrability 
assumptions in and +oo (namely /j^max ^j^, e'^'^^ ^( do;) < oo for some constant 

oo;. 

If we choose h such that h-g-C with g{-,oj) - q{-,ijj)f{-,uj), we have the following: 

• By construction of /, using ()3.18p in (13.170 . T{h, I) - \\l ||^ -i q = 1) 

• II ^ 11^ 2 ^ ^ Is f'^^ ^ ^- ^^P|| ?»|| 2 1^1 1^(^' 01 ^ ^) where C is indepen- 
dent of I e Q such that |U ||„ _i „ = 1- 

Applying Proposition 13. 7( we obtain the existence of some pi e n O such that 
r(pi, ~ /(Oi I & G- But by density, this is also true for I e H^^^ n C □ 
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Proof of Lemma \3.(k Since there exists some constants C, c > such that for all oo e Supp(;u), 
^ e S, < c ^ q{9, u) ^ C, £ is continuous on H^^^nO (as well as on L^nO). More precisely, by 
Riesz theorem, there exists a unique e e H^^^nO such that for all I e H'^^nO, - (e , I) ^ _^ ^. 
One can be more explicit: a simple calculation shows that this {0,u) e{6,uj) corresponds 
to the primitive £{9,u}) - -q{9,u}) cos{6), that is: 

V6I e S, w e Supp(/i), e(e, oj) ^ -deqie,uj) cos{e) + q{e,uj)sm{e) . (3.19) 

Let us introduce the following function p2 e n O: 

-B(e,Lo) f~ f~g 

P2(^,c^) . / eS(«.-')+4-'^ dn + / e-^("''^^-^(^'"'^ dtx, (3.20) 

1 6 ^ s *y 

for 

5(6*, uj) = -2 {Kr {cos(e) - 1) + ^^6*) . 
Then one readily verifies that Lp2 is proportional to e. □ 

4. Global spectral properties of operator L (Proof of Theorem 12. 7|) 

The purpose of this section is to prove Theorem 12.71 The main idea of the proof is to 
decompose the operator L defined by (j2.ip on the domain T> given by (|2.2p into the sum 
of a self-adjoint operator A (in a weighted Sobolev space for appropriate weights) and a 
perturbation B which will be considered to be small w.r.t. A. Namely, it is straightforward to 
verify from (|2.ip that L can be written as L - A + B where, for all h €D, for all uj e Supp(/i), 

Ah{e,u;) ^d^eHO,^) - de(h{e,u){J * qo){d) + qo{e){J * (4.1) 

and, 

Bh{e,u;) -de[h{e,uj){{J * {q - qo))^.{e) + uj} + {q{e,u) - qo{e)){J * %). (4.2) 

We divide the proof of Theorem 12.71 into three parts: in § 14. H we prove that A is essentially 
self-adjoint (and thus generates an analytic semigroup) in some weighted-Sobolev space (recall 
§ 12. 2p for an appropriate choice of weights. Note that this section strongly relies on the fact 
that /i is a binary distribution (recall (12. 4p ). 

The purpose of § 14.21 is to establish precise control of the size of the perturbation B w.r.t. 
A. The third step of the proof, carried out in § 14. 3^ consists in deriving similar spectral 
properties for the operator L = A + B, especially the fact that the spectrum of L lies in the 
part of the complex plane with negative real part. 

4.1. Spectral properties of the operator A. In this paragraph, we prove mainly that A 
defined in (j4.ip is essentially self-adjoint for a Sobolev norm that is equivalent to the norm 
II -IIh^ defined in S [2X2l 

Since we are working in the domain P (recall (12. 2p ). the test functions h are such that 
j-^h{-,oj) dji = I (/i(-,+a;o) +/i(-,-a;o)) has zero mean value on S. The idea of this paragraph 
is to reformulate the operator A in terms of the sum |(/i(-, +ojq) + h{-, -wq)) and the difference 
i(/i(-, +ujo) ~ -ujo)); namely, we define the following 2x2 invertible matrix: 
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and for /i e P, let ( " ) := M • /i, namely 

I n(-) := i(/i(-,+cJo) + /i(-,-Wo)), 
1 := i(/l(-,+tJo) -/l(-,-Wo))- 

We are now able to define the following operator: A:- M oAo M^^ , defined on the domain V 



(4.3) 



V:^\{u,v)eC^{S)xC' 



de^o 



given by 



V{u,v) e V, 




^ ^dQ2U - de [n( J * qq) + qo{J * u)]^ 
1 



^ -de^v - de[v{J * qq)] 



(4.4) 



(4.5) 



/ 



The remarkable observation is that operator A is now uncoupled w.r.t. variables u and v] 
consequently, in order to diagonalize A^ it suffices to diagonalize both components of A, 
namely Ai and ^2- This is the purpose of Propositions 14. l l and 14.21 

We make use here of the weighted Sobolev norms || • defined in (j2.7p for different choices 
of k{-): an important remark is that the first component Ai (with domain {u e C^(S), J^u - 
0}) is exactly equal to the McKean-Vlasov operator Lg^ with no disorder defined in (j2.1ip . 
Following § 12.31 and [5] , the natural space in where to study the operator Ai is the space 
U-^ defined in ([22|, for the weight fc(-) = qo{-) the non-disordered stationary solution (recall 
(ll.lSp ). In this space, we have 

Proposition 4.1. In H~^^, Ai is essentially self-adjoint with compact resolvent and spectrum 
in the negative part of the real axis. is a one- dimensional eigenvalue, spanned by dgq^. The 
spectral gap Xk{Ai) - XxiLq^) between and the rest of the spectrum is strictly positive. 

Moreover, the self-adjoint extension of Ai is the infinitesimal generator of a strongly con- 
tinuous semi-group of contractions T'i(t) on . For every < a < ^, this semigroup can be 
extended to an analytic semigroup Ti{z) defined on Aq, = {A; |arg(A)| < a} and one has the 
following estimate on its resolvent (where Sq, = {Ae C; |arg(A)| < f + a}u{0};.- 



Va e 



(0,f), VA 



eS,, \\RiX,Ai)\ 



-1,90 



1 



1 



l-sin(a) |A| 



(4.6) 



The second component A2 is a second order ordinary differential operator, with domain 
C^(S). As we will see in § 14.1.21 the natural space in which to study A2 is the space H~^, for 
the choice of the weight function 9 w(9) - 



„-*(e) 



with 



^{9):^-2Krocosi9), (4.7) 
where tq is given by ()1.19p . Namely, we have 

Proposition 4.2. The operator (^25^^(8)) is essentially self-adjoint in H^^ and has compact 
resolvent. Hence, its spectrum consists of isolated eigenvalues with finite multiplicities. The 

kernel of A2 is of dimension 1, spanned by w(9) = 7 — —. Moreover, we have the following 

Js 

spectral gap estimation: 



V^;eC2(S), -[A2V,v)_^^^ 



-4Kro 



W 



-1,10 



(4i 
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SO that the spectrum of A2 lies in the negative part of the real axis and the distance between 
and the rest of the spectrum Xk{A2) is at least equal to ^—^ — ■ One also has explicit estimate 
on the resolvent of A2: 

Va6(0,^), VAeS,, || i?(A, ia) ||_, ^ ^ (4-9) 

2 l-sin(a) |A| 



Putting things together, the natural norm for the operator A - (Ai,A2) is the Hilbert- 

1 ^ 

norm, for {u,v) e V: {\\u\\-iqg + II ^ ll-i «;) ^ • But since A is the conjugate of A through the 
invertible matrix M, to say that A is essentially self-adjoint for the previous norm is equivalent 
to say that A is essentially self-adjoint for the corresponding conjugate norm: 



V/i €T>, II /ill „ 



^ {h{-,+ujo) + h{-,-ujo)) 



2 

+ 

-1,90 



^ {h{-,+u}o) - h{-,-ujo)) 



2 



-l,ui J 



(4.10) 

The results of § 14.11 can be summed-up in the following proposition, which is an easy conse- 
quence of Propositions 14.11 and 14. 2t 



/a.(0,|) 



Proposition 4.3. For the norm ||-||j^^ defined in (j4.10p . the operator (A,T>) is essentially 
self-adjoint, with compact resolvent. In particular, the spectrum of (the self-adjoint extension 
of) A is pure-point, and consists of eigenvalues with finite multiplicities. Moreover it lies in 
the negative part of the real-axis and A is the infinitesimal generator of an analytic semigroup 
of operators Ta{z) defined on a domain - {z € C; \ arg(2:)| < a}, for any < a < ^. One 
also has the following estimate about the resolvent of A: 

VAeS,, ||i?(A,^)||^^ ^ (4.11) 

l-sm(a) |A| 

The kernel of A is of dimension 2, spanned by Ideqo + / -^^ idelo ~ 1 CLnd the eigenvalue 
is separated from the rest of the spectrum with a distance \k{A) ■■- min (A/^ ( Jii), Ai<-( ^12)), 
where XxiAi) and Xk{A2) are defined in Propositions \4-l\ and \4-2\ respectively. 

Remark 4.4. Note that this norm \\ ■ is equivalent to the norm \\ ■ ||jj^^ defined in S \2.2.^ 
since the weights qo and w are bounded above and below. In particular, in the space Hj,^ , the 
operator A (although no longer self-adjoint) still generates an analytic semi-group with the 
same spectrum and the same spectral gap. 

The aim of paragraphs § 14.1.11 (resp. § I4.1.2p is to prove Proposition 14.11 (resp. Proposi- 
tion 1121). 



4.1.1. Spectral properties of A\: proof of Proposition Much about the operator Ai has 
already been done since it exactly corresponds to the linear evolution operator L^g of the 
non-disordered Kuramoto model studied in [5j: we know from Proposition 12.41 that Ai is 
essentially self-adjoint and dissipative in Hg^. It remains to prove that Ai generates an 
analytic semigroup Ti{t) in an appropriate domain. We refer to classical references [20 11271133] 
for detailed definitions of analytic semigroups of operators defined on a sector of the complex 
plane. In particular, we recall the result we use concerning analytic extensions of strongly 
continuous semigroups. Its proof can be found in [33l Th 5.2, p. 61]. 
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Proposition 4.5. Let T(t) a uniformly bounded strongly continuous semigroup, whose infin- 
itesimal generator F is such that 6 p^F) and let a e (0, -1). The following statements are 
equivalent: 

(1) T{t) can he extended to an analytic semigroup in the sector A.^ - {A e C; | arg(A)| < a} 
and \T{z)\ is uniformly hounded in every closed suh-sector Aa', ol < a, of Aa, 

(2) There exists M > such that 

p{F) D = |a e C; I arg(A)| < ^ + «} u {0}, (4.12) 

and 

M 

\\R(X,F)\\^—, AeE,A*0. (4.13) 



We are now in position to prove the rest of Proposition I4.lt we know from [5^ Prop. 2.3, 
Prop. 2.6] that for any A > 0, X - Ai is positive with range H^^ . Consequently, we are in 
position to apply Lumer-Phillips Theorem (see |33[ Th 4.3 p. 14]): Ai is the infinitesimal 
generator of a Cq semi-group of contractions denoted by Ti(t). 

The rest of the proof is devoted to show the existence of an analytic extension of this 
semigroup in a proper sector. We follow here the lines of the proof of [33\ Th 5.2, p. 61-62], 
but with explicit estimates on the resolvent: let us first replace the operator Ai by a small 
perturbation: for all e > 0, let Ai^^ := Ai - e, so that belongs to p{Ai^s). The operator 
Ai^£ has the following properties: as Ai, it is self-adjoint and generates a strongly continuous 
semigroup of operators (which is Ti^e(t) = Ti{t)e^^^). Moreover, since Ai^^ is self-adjoint, we 
have 

VAeCxR,||i?(A,ii,,)||_^^^^ ^^-^^^ 
and since the spectrum of Ai^f, is negative, for every A e C such that 1H(A) > 0, we have 

(4.15) 



Let us prove that for A e 



RiX, ii,e) ||_, ^ ■ ^ . (4.16) 

i-,go l-sm(a) |A| 



Note that (j4.16p is clear from (j4.14p and (j4.15p when A is such that 9\{X) ^ 0. Let us prove it 
for A e Sq, with 1K(A) < 0. Consider o" > 0, r e R to be chosen appropriately later and write the 
following Taylor expansion for i?(A,j4e) around a + ir (at least well defined in a neighborhood 
of a + ir since o" > 0) : 

oo 

i?(A,ii,,) = ^ii(a + zT,ii,,)"^H(^^ + ^r)-A)". (4.17) 

n=0 

This series R{-,Ai^^) is well defined in A e with 9^(A) < if one can choose a, r and 
k e (0, 1) such that || R{a + iT,Ai^^) || _^ |A - (fi + iT)\ ^ k < 1. In particular, using (j4.14p . it 

suffices to have |A - (a + zt)| ^ /c|t| and since cr > is arbitrary, it suffices to find k e (0, 1) 
and T with |A - «t| ^ /c|t| to obtain the convergence of ()4.17p . For this A e Eq, with 9^(A) < 0, 
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let US define A' and r as in Figure HI Then, |A - ir| ^ |A' - ir| = sin(a)|T| with sin(a) e (0, 1). 
So the series converges for A e and one has, using again (j4.14p . 

||i?(A,ii,,)||_ ^ / ^ (4.18) 

i-,go (1 - sin(a))|r| 1 - sin(a) |A| 

The fact that Ti^i,(t) can be extended to an analytic semigroup T-i^^(^z^ on the domain 




Figure 4. The set Sq. 

is a simple application of (j4.16p and Proposition 14.51 with M . Let us then define 

Ti(z) ■- e^^Ti^i;(z), for z e so that Ti(z) is an analytic extension of Ti(t) (an argument of 
analyticity shows that Ti(z) does not depend on e). 

Note that estimation (|4.6p can be obtained by letting e in ()4.16p . □ 

4.1.2. Spectral properties of A2: proof of Proposition ^2 may be written as 

A2V^^dlv + de{vKrosm{-)), (4.19) 

where ro = ^'o(2irro) (recah ([119]) ). 

One recognizes in JL2 a second order differential operator on C^(S) which can be identified 
to a simple instance of a Fokker-Planck operator with a sine potential (see [35] )• This operator 
can easily be seen, by integrations by parts in an appropriate weighted L^-space as a Sturm- 
Liouville operator acting on C^, 27r-periodic functions. A vast literature exists on the spectral 
properties of such operators (see [SHl ESI El EU O [M] ) . But the problem is that a L^-norm is 
not appropriate for the future study of the SPDE (II. 7p : a look at the covariance structure of 
the noise W (see p.9p ) shows that W naturally lives in a H^^-space instead of a L^-space. 

An easy calculation shows that A2 can be rewritten in terms of the weight function $ 
defined in (lO) : 

A2V^^de{e~''de{e''v)). (4.20) 
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Let w be: 



is ^ 

One directly sees from (j4.20p that w lies in the kernel of A2- A2W = 0. 

We are now in position to prove Proposition 14.21 let us place ourselves in the framework 
of the weighted-Sobolev spaces (L^,(-, •)2^„) and (H^^,(-, constructed in S 12.2.1^ in 

the particular case of k{-) - w{-). 



Proof of Proposition \4-3\ In H^^, the operator (j42,C^(S)) is formally symmetric: indeed, for 
n, u e C^(S), for tio ad vq defined by (12.6j) . we have successively, 



^ e^uvQ, 



r!m>,_i r r 1 rum, ^ 

2 Js w 2 Js Js 2 Js w ^ ' 
The latter expression being symmetric in u and v. 

Let us prove that (Ji2,C^(S)) is essentially self-adjoint: let S2 be the following Dirichlet form 

£2{n,v):^{u,{l-A2)v) ^ ^ [ufv^ f^^lf^. (4.23) 

^ ^ ^ ^ ' '-I''" Js Js Js w 2Js w ^ ^ 

Then it is easy to see that £2 is a continuous bilinear form on (thanks to Poincare 
inequality). Moreover £2 is coercive, namely, for all u e L^: 



^^IkllL- (4.24) 



Since for all / e H^^, the linear form v {v , /)_^ ^ is continuous on L^, an application of 
Lax-Milgram Theorem shows that for such an / e H^^ there exists an unique u e such 
that for all v e 

£2{v,u)^{v,f)_,^^. (4.25) 
It is then easy to see that /g / = fgU and that for almost every 6 € S, 



r^.r^. (4.26) 
Jo w Jo w 



2 w{e) 

Since u e L^, Z//o admits a C^-version and if we assume that / is square-integrable, the same 
argument holds for the first term of the right-hand side of ()4.26p . So, if / is square integrable, 
no admits a C^-version. To sum-up, if we suppose that / is continuous, there exists u e C^(S) 
such that, applying dg (e"*96»(0) to (|4.26p : 

/ = /o + / j ^« = ~de {e'^de (e*no)) + uo+Pw{u)w, 

^-A2Uo + u^{l-A2)u. (4.27) 

But since those functions / are dense in , we see that the range of 1 - is dense so that 
A2 is essentially self-adjoint. 



26 ERIC LUgONi 

Secondly, the spectral gap estimation (14. Sp holds: indeed, for every u e C^(S), we have 
using (j4.22p and Poincare inequality: 



$ 2 



2 \Js I 



-l,ui 

Moreover, A2 has compact resolvent: it suffices to prove that X- A2 has compact resolvent 
for at least one value of A. We prove it for A = 1 which is indeed in the resolvent set, 
thanks to the beginning of this proof. For u e H^^, let us consider / := (1 - A2)~^u so that 
(/, (l-^2)/)_^^ = (/, "u)-!^- Using the coerciveness of £2, one has for some constant c, 

c||/ll2,«, ^ (/' ^ II / II • Using the continuous injection of into H^^ 

(say II • ||_i ^ C II • II2 some positive constant C), one has 

II/II2,. ^fll^lLi,^,- (4.28) 

So (1-^2)^ maps sequences that are bounded in H^^ into sequences that are bounded 
in L^. It remains then to prove that the injection of into H^^ is compact. This 
is indeed true since for every v e H^^, one has, by Cauchy-Schwartz inequality |Vo(6') - 
Void')] ^ C II Il2,«; V|6'-6"| ^ C II w II 2^^ \/\0-0'\. That means that, by Ascoli-Arzela The- 
orem that the sets [v e H^^; || f II2 t„ ^ est} are relatively compact in C(S) and also in L^. 
That completes the proof. 

The fact that A2 generates an analytic semigroup 22(2:) on the same sector as well as 
estimation ()4.9p can be derived in the same way as in § 14.1.11 That concludes the proof of 
Proposition 14.21 □ 



4.2. Control on the perturbation B. In order to derive spectral properties for the operator 
L - A + B, we need to have a precise estimation about the smallness of the perturbation B 
w.r.t. operator A studied in the previous paragraph § 14.11 

Remark 4.6. Note that for simplicity, we work now with the norm || • ||fj (recall Remark l2.3p : 
as already mentioned this norm is equivalent to the norm || • used in § 14.11 Recall also 
the definition of the space (L^, || • Hg) defined in Remark 12.21 and of (L^j || ■ II^ 2) defined in 
Remark 12.31 

Secondly, since the whole operator L is no longer symmetric in H^^, its spectrum need not 
be real. Thus, we will assume for the rest of this document that we work with the complexified 
versions of the scalar products defined previously in this paper. The results concerning the 
operator A are obviously still valid. 

The smallness of the perturbation B with respect to A can be quantified in terms of the 
difference || q{-.,uj) - (7o(") Hoc ^ Supp(/i). For the ease of exposition, we do not attempt to 
derive precise estimations of this difference || q{-.,oj) - goC") ||oo (Lemma 14. 7p and of coefficients 
a{oJo) and 6(wo) (Lemma 14. 8p . in terms of the coupling strength K. We will denote as c a 
positive constant (depending on K) which may change from a line to another. 
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Lemma 4.7. For bj> {) and K > 1, let us define 

\\q-qo\\oo--- sup \qie,u)-qoie)\. (4.29) 

6»eS,|n| ^cj 

Then || ^7 - (Zo II cxj ~ 0{uj), as u ^ 0. 

Proof. This is clear since one can bound d^q{9^uj) uniformly in (0,a;), as a; ^ (by a constant 
depending on K). □ 

The purpose of the following proposition is to prove the boundedness w.r.t. A of the 
perturbation B in terms of the size of the disorder: 

Proposition 4.8. The operator B is A-bounded in the sense that there exist positive constants 
a{uJo) - a(a;o,i^) and ^(wo) = b{(jjQ,K) such that 

V/i e V, II Bh IIh^ ^ a(wo) II h \\^^ + b{uo) \\ Ah \\^^ , (4.30) 

and moreover, for fixed K > 1, 

a{tJo) - O^ujq), and b{ijJo) - 0{ijJo), asojQ-^O, (4-31) 

Proof of Proposition ]^-.^ Recall that in what follows, {h)^(-) - /i(-, a;)/x( dw) is the aver- 
aging of h(-,u) and that c is a constant (depending on K) that may change from a line to 
another. The proof consists in two steps: let us prove firstly that there exist some constant 
OiK,u}o such that for all h €V, 

\\Bh\\^^^aK,ojo\\h\\f,^2^ (4.32) 

Indeed, for given h € V, for all u e Supp(/z), we have ||i?/i(-, (j;)||_i = || i3/i(-, w) Hg, where 
Bh{-,io) is the appropriate primitive of Bh{-,uj) in H~^, (recall Remark I2.2p namely: 

Bh:= -h{{J *{q- qo))f, + u})-{q- qo) ■ (J * h)^ 

Using the boundedness of go and the bounds | (J * e) | ^ 4K || e and |(J*/i)^| ^ ^ || (/i)^ Hg, 
it is easy to deduce that, for some constant c> 0: 

m ^c{\\q-qo L +a;o)(|/i| + || (% II2). (4.33) 

Consequently, 

II Bh IIh^ = ((II Bh ^ c(|| g - go L + ^0) II h \\^^^ , (4.34) 

so that (|4.32p is satisfied for some coefficient otK,u)(, such that (Lemma 14. 7p oix,u]o ~ 0(^ujq) as 
Wo 0. 

The second step of the proof is to control the L^-norm || /i ||^ 2 °^ ^ with the H~^-norms of 
Ah and h: namely we prove that there exist constants '^k and 5k such that 

II h II ^ II Ah II +5K\\h II . (4.35) 

The proof is based on a usual interpolation argument: for all integer n > 1, for any / 6 C^(S), 
one has 



^Vn / 2 ^^77 / 2 + ^^ • (4-36) 

\/n 1 In 
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Let us use this interpolation relation (|4.36p to derive (I4.35P : for all h € V, u € Supp(/i), one 
has 

II 2 



^(•>'^) II2 = [^^(•'^)] +\\ho{-,u})\ 
Applying relation ()4.36p with /(•) = ^o('i^) we obtain 



h{-,u)) 



n 



deh{-,uj) II 2 
2n 



where we used the fact that dgho{-,iLi) - dgh{-,u). Integrating w.r.t. /x, 



h\ 



M>2 



f f 

Jr Js 



d/i + 



n 



II 2 

\\fj.,2 



no 11: + 



deh\ 



fj.,2 



2n 



(4.37) 



(4.38) 



(4.39) 



But, as previously for the operator B, a simple calculation shows that for all u e Supp(/x), we 
have ||A/i(-,(x>)||_i = II Ah{-,Lo) II2, where Ah is the appropriate primitive of Ah in H""*^ (recall 

my- 

(4.40) 
(4.41) 

(4.42) 



Ah = ^dgh - h{J * qo) - qo{J *^)^^ + (H^ * Qo) + qo{J * h)^) 
so that, for some constant c > 

\\9eh\\l^2 ^ 12||^/i||J^2 + c||/i||^_2- 
Injecting this inequality in ()4.39p . one obtains, for some constant c > 



11^,2 



^ x u ^ i " + ^ (12 II Ah 11^,2 +c\\h 



2 

Choosing n > 1 sufficiently large so that the coefficient in front of || /i || ^ 2 right-hand 
side of ()4.42p is lower than ^ leads to (for some constant c> 0): 

II ^ £2 ^ 2 ^ /i d^ + c II £,2 + c\\Ah ||^_2 ^c\\h f^^ + c\\Ah ||^^ , 

which shows (j4.35p . Putting together estimates (j4.32p and (j4.35p . we find the A-boundedness 
of B (j4.30p with coefficients a(cjo) and 6(a;o) which satisfy (j4.3ip . thanks to Lemma 14.71 □ 

Proposition 4.9. The operator B is A-compact, in the sense that for any sequence {hp)p ^0 ^ 
T>^ such that \\ hp ||jj and \\ Ahp ||jj are bounded, there exists a convergent subsequence for 
(Bhp) 



Proof of Proposition \4-9[ Let {hp)p ^ a sequence in V such that || hp ||jj^ and || Ahp ||jj^^ are 
bounded by a constant c. A closer look at the operator B defined in (j4.2p and the definition 



of the norm 



in (|2.8p shows that it suffices to prove that there exists a subsequence 



p ii/i 2 ^ + c II hp 11^ 2) so that 



(^Pfc)fc^o that is convergent in L^. In particular, for all p ^ 0, || ^/ip ||jj ^ c. Using this 
boundedness and (j4.4ip . we have || dgk 

\\dehpl^,^c.c(f^\f^ 

< c + c 



dfi + II ho,p 



Jr fs 



n 2 



dfjh. 



2n 



(4.43) 
(4.44) 
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where we used again (j4.36p for / = T-Lqp^-jUj). Choosing a sufficiently large n > 1 leads to 
II dgho^p 11^ 2 ~ II ^^^p II2 ^ '^'^'^ II II ^ ^ constant c independent of p ^ 0. An easy appli- 
cation of Cauchy-Schwartz inequality leads to |/io,p(^, w) - /io,p(^') '^)| ^ II dgho^p ||^ 2 \/|^ - 
for all cj e {±^0}- Since the functions {hQ^p)p ^ are such that /g /io,p = for all 0, Ascoli- 
Arzela Theorem and the previous bound show the existence of a convergent subsequence 
(/io,pfc) (for each uj e Supp(;u)) in the space of continuous functions on S. In particular, this 
subsequence is convergent in and is renamed (/io,p)p^O) with a slight abuse of notations. 

The fact that \ /g hp\ dfj, ^ c shows that one can extract a further subsequence of {hp)p j, 
which is also convergent in L^. This concludes the proof. □ 

4.3. Spectral properties of L - A + B. We are now in position to extent the spectral 
properties we obtained in § l4.1l for the operator A to the perturbed operator L - A+B, relying 
on the control we have the perturbation B (see § 14. 2p and on the theory of perturbation of 
operators (p3j) and analytic semi-groups ([33]). Recall again that the disorder is supposed 
to be binary (see (j2.4p ). 

4.3.1. The spectrum of L is pure point. 
Proposition 4.10. For all K, for all ujq > 0, 

(1) the operator (L,T)) is closable. In that case, its closure has the same domain as the 
closure of A, 

(2) the closure of L has compact resolvent. In particular, its spectrum is pure point. 

Remark 4.11. Note that Proposition \4-10 is valid without any assumption on the smallness 
of loq, since it relies on the relative compactness of B with respect to A (Prop. \4-9^ - 

Proof of Proposition \4^7d It is a simple consequence of the relative compactness of B w.r.t. 
the self-adjoint operator A. The first assertion of Proposition |4TT0] is a consequence of |23l Th. 
1.11 p. 194] and the second assertion can be found in [28^ Lemma 3.6, p. 17] for example. □ 

4.3.2. L generates an analytic operator. We prove mainly that the perturbed operator L 
generates an analytic semigroup of operators on a appropriate sector. An immediate corollary 
is the position of the spectrum in a cone whose vertex is zero. We know (Proposition 14. 3p 
that for all < a < ^, ^ generates a semigroup of operators on - {z € C; \ arg(z)| < a}. 

Proposition 4.12. For all K > 1, e > and < a < ^, there exists wi > (depending on 
a, K and e) such that for all < wq < wi, the spectrum of L lies within the sector Qe,a 
{a e C; ^ + a ^ arg(A) ^ ^ - a} u {A e C; |A| ^ e} . For such ujq, L generates an analytic 
semigroup of operators on Aq,', for some a' e (0, ^). 

Proof of Proposition \4^T^ Let < a < | be fixed. Thanks to (|4.1ip . there exists a constant 
c > (which comes from the equivalence of the norms || • || jj^ and || • || ji/^ ) such that for every 

A e := {a e C; I arg(A)| < f + a}, we have the estimates: 

II^(^'^)IIh. ^7^ ^^^^ ^^d, \\AR{\,A)\\^^l + - . 

(l-sm(a))|A| (l-sm(a)) 

Then for A e S^, /i e 2?: 

II BR(\, A)h IIh^ ^ a(a;o) II ^(A, A)h \\^^ + b{uo) \\ AR{X, A)h \\^^ , 

\ (l-sm(a))|A| \ l-sm(a)// 
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Let us fix e > and choose uji such that: 



max 



( 4K..)(l^T^))<l. (4.45) 

For this choice of wi, for ah < c^o < wi, for any A e Sq, such that |A| ^ e ^ i~s\n{a) ' 
have II -Bi2( A, 74)/i II jj ^ III^IIh ' ^"^^ thus the operator \ - BR{\,A) is invertible with 
II 1 - BR(X, A) II ^ 2. Since it can easily be shown that 

R{X,A + B)^R{X,A){l-BR{X,A)y\ 
one deduces the following estimates about the resolvent of the perturbed operator L - A + B: 

VAeS,, |A|^e, || ii(A, L) ||h^ ^- (4.46) 

(l-sm(a))|A| 

The fact that the spectrum of L lies within Qe^a is a straightforward consequence of (I4.46p . 
Secondly, (I4.46P entails that L generates an analytic semigroup of operators on an appropriate 
sector. Indeed, if one denotes by := L-e, one deduces from (I4.46P that e p{L2e) and that 
for all A 6 C with 9^(A) > (in particular, |A| < |A + 2e|) 

ll^(A,i^2.)||„ = ||i?(A + 2e,L)||„ ^ — 1— — , 

^ (1 - sm(a))|A + 2e| 

2 

^ (l-sin(a))|A| ' ^^'^"^^ 

Hence, using the same arguments of Taylor expansion as in the proof of Proposition 14.11 and 
applying Proposition 14. 5| one easily sees that generates an analytic semigroup in a (a 
priori) smaller sector Aq,', where a' e (0, ^) can be chosen as a' ■- | arctan ^ ^— ) • But if 
L2e generates an analytic semigroup, so does L. □ 

4.3.3. is an eigenvalue with multiplicity 2. Let us fix > 1, a e (0, p e (0, 1) and define 
£ - pXk{A). Applying Proposition 14.121 we know that for small loq (depending on K, a, p), 
L generates an analytic semigroup on Qe^a- Let 0^^ := {A e @e,a', ^{^) ^ 0} be the subset of 
Qs,a which lies in the positive part of the complex plane (see Figure [5|). 

The purpose of this paragraph is to show that one can choose a perturbation B small 
enough so that no non-zero eigenvalue of A + B remains in the small set @^ a- 

To do so, we proceed by an argument of local perturbation: we know that the distance 
Xk{A) - min (A/f ( Ai), Ai<-(^2)) between the eigenvalue and the rest of the spectrum of A 
is strictly positive. In particular, one can separate from the rest of the spectrum of A by a 
circle centered in with radius (^)Aa'(^)- Note that the appropriate choice of e made 
at the beginning of this paragraph ensures that the interior of ^ contains Q^^a (Figure [5]). 

The main argument is the following: by construction of is the only eigenvalue (with 
multiplicity 2) of the non-perturbed operator A lying in the interior of A principle of 
continuity of eigenvalues shows that, while adding a small enough perturbation B to A, the 
interior of still contains either one eigenvalue with algebraic multiplicity 2 or two eigenvalues 
with multiplicity 1; those perturbed eigenvalues remain close but are a priori not equal to 
the initial eigenvalue (see Figure [5]) . 

But we already know that for the perturbed operator L - A + B, is always an eigenvalue 
(since Ldeq - and Lp - dgq, recall Th. 12. 5p . Therefore, the algebraic multiplicity of for 
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(a) Position of the spectrum for (b) Possible position of the spec- 

the self-adjoint operator A trum for the operator L 



Figure 5. The domains Qe,a and G^^ (in light grey). Note that the two 
dimensional eigenvalue for the operator A (I5ap may split in two single eigen- 
values for the perturbed operator L ()5bp . These eigenvalues are the only ones 
within the circle 'if. But since we already know that the eigenspace in of L 
is of dimension 2, is still a double eigenvalue for L, by uniqueness. 



the operator L is at least 2. By uniqueness, one can conclude that is the only element of 
the spectrum of L within and is an eigenvalue with algebraic multiplicity exactly 2. In 
particular, there is no element of the spectrum in the positive part of the complex plane. 

In order to make this argument precise, we need to quantify the appropriate size of the 
perturbation B, by explicit estimates on the resolvent R{X,A) on the circle 'if: 

Lemma 4.13. There exists some explicit constant o^{K) only dependent on K, such that 
for all X^'^ , 

\\R{\,A)U^^c^,{K), (4.48) 
II AR{\A) IIh^ 1 + .c^^K). (4.49) 

One can choose ccg{K) as c^g{K) - xJ(a) ( p+T ' ) • 

Proof of Levama \4-l'J\ Applying the spectral theorem (see |12j Th. 3, p. 1192]) to the essen- 
tially self-adjoint operator A, there exists a spectral measure E vanishing on the comple- 
mentary of the spectrum of A such that A - f-^XdE(X). In that extent, one has for any 



(4.50) 



32 



ERIC LUgONi 



In particular, for C e ^ 

1 max I -jy , -r^ ) 
\\RiC,A)\\^^^ sup <^flL±Rl, (4.51) 

The estimation (j4.49p is straightforward. □ 

We are now in position to apply our argument of local continuity of eigenvalues: follow- 
ing [231 Th III-6.17, p. 178], there exists a decomposition of the operator A according to 
H^^ - Fq® 1^0 (in the sense that ^4^0 c Fq, AF^q c i^g and PqV c V, where Pq is the projec- 
tion on Fq along F_q ) in such a way that A restricted to Fq has spectrum {0} and A restricted 
to i^o has spectrum cr{A) \ {0} £ {A e C; 9^(A) < 0}. Let us note that the dimension of Fq 
is exactly 2, since the characteristic space of A for the eigenvalue is reduced to its kernel 
which is of dimension 2 (see Prop. 1^3]) . 

Then, applying [23, Th. IV-3.18, p. 214], and using Proposition 14.81 we find that if one 
chooses UJ2 > 0, such that 

sup (a(w2) II R{X,A) ||„ + 6(w2) II AR{X,A) ||„ ) < 1, (4.52) 

then for all < (x>o < i^2) the perturbed operator L is likewise decomposed according to 
H^^ = Go ® G<o, in such a way that dim(Fo) = dim(Go) = 2, and that the spectrum of L is 
again separated in two parts by . But thanks to Theorem 12. 5^ we already know that the 
characteristic space of the perturbed operator L according to the eigenvalue is at least of 
dimension 2 (since Ldgq - and Lp - dgq). We can conclude, that for such an < (x>o < i^2) 
is the only eigenvalue in ^ and that dim(Go) is exactly 2. 

Applying Lemma 14.131 we see that condition (j4.52p is satisfied if we choose u;2 > so that: 

a{uj2MK) + 6(^2) (1 + (^) o^iK)^ < 1. (4.53) 

In particular, in that case, the spectrum of L is contained in {z e C; yi{z) ^ 0}. 
Theorem 12.71 is proved. 



5. non self-averaging phenomenon for the fluctuation process (proof of 

Theorem \T9\) 

The purpose of this section is to prove our main result Theorem 12.91 that is the linear 
asymptotics (I2l9]) for the SPDE (fT7fl> . 

In our framework, (recall that /i = |(^-wo + <Jwo))i *he solution rj of evolution (|1.7p in 
S'(S X R) acts on test functions ip of the form ip - ((p(-,+ujQ),(p(-,-uJo)). In particular, one 
can understand rj as an element of H^^ by identifying rj with (T/(^Q,?7-ajo), where, for any 
smooth function • S R, r]i^g{ip) 7]{ip,0) and r]-u)o{ip) '■= ??(O,'0). Defining analogously 
14^±wo the Wiener process W in ()1.9p . the object of interest is then 



yt>0, rit^m+ f Lri ds + Wt- 
Jo 



(5.1) 
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5.1. The noise 1^ as a cylindrical Brownian Motion. We first focus on the regularity 
properties of the noise W: in the stationary case {qt - qt\t=o - Q for all t > 0) the covariance 
defined in (|1.9p becomes, for any regular functions on S x R, ipi and ip2, s,t> 0: 

BiWsi^i)Wtiip2))-isAt) f f deipideip2qdXdf,. (5.2) 

Jr Js 

Consequently, it is easily seen that (W+o^o,*) W^-wo,t) is a couple of two independent Gaussian 
processes with covariance (where ipi,ip2 ^ S -f R) : 

Vw6{+wo,-^o}, ^{W^A^iWuM^)) - \{s At) J^deAdei^2q{-,u:). (5.3) 

In what follows, we will denote by Hq the closed subspace of H^^ consisting of elements of 
H^^ with zero mean- value; in particular the norm ||-||fj defined in ()2.8p coincides on Hq 
with: 

V/^eHo, ll^llH,-(2vrJ^^^')', (5.4) 

where we recall that T-L is the primitive of h such that f^Ti - 0. The purpose of this paragraph 
is to show that, following [9l p. 96], W has the same law as a Q- Wiener process in the Hilbert 
space Ho, for an appropriate bounded symmetric operator Q on Hq: indeed, if one denotes 
by X a Q- Wiener process on Hq, with the following definition of Q 

V/ieHo, Qh-.^deiqn), (5.5) 

then one readily verifies that the Gaussian process iW(ip))ip has the same law as (X(ip))^ 

Note that if the weight q had not been present in (15. 3p . it would have been straightforward 
to see that is a Q- Wiener process on Hq, with Q - I. The fact that this weight is present 
entails some technical complications, but one has to consider the operator Q defined in (15. 5p 
only as a perturbation of the case Q - I. 

5.2. Existence and uniqueness of a weak solution to the fluctuation equation. We 

now turn to the existence and uniqueness of a weak solution of (j5.ip . We recall that any 
H~^-valued predictable process 7]t, t e [0,r] is a weak solution of (|5.ip if the trajectories of r] 
are almost surely Bochner integrable and if for all if e T){L*) and for all t e [0,r] 

miv) ^ mif) + £ VsiL*ip)ds + Wtiip). (5.6) 

Proposition 5.1. Equation has a unique weak solution in H~^, given by the mild formulation 

m-TL{t)rio + J^TL{t-s)dWs, te[0,r]. (5.7) 

To prove Proposition [5Tl one needs to define properly the stochastic convolution WL{t) := 
/o TL{t - s) dWs- In this purpose, let use prove firstly that the inverse of A is of class trace: 

Lemma 5.2. The operator A^^ is of class trace in H^-*^. Equivalently, if {xi^^)n^ i is the 
sequence of eigenvalues of the self-adjoint operator A, one has 

oo -I 
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Proof of Lemma \5.S[ Since A and A - M o A o M^^ (recall ()4.5p ) are conjugate, it suffices to 
prove (|5.8p when A is replaced by A. The idea of the proof is that identity ()5.8p is true when 
A - (Ai,A2) is replaced by (-A,-A) and that A is only a relatively-bounded perturbation 
of this case. More precisely, the proof relies on the following MinMax Principle [13, p. 1543]: 

Proposition 5.3. Let {F,T>{F)) a self-adjoint linear operator on a separable Hilhert space % 
such that F is positive, with compact resolvent. We denote by S"' the family of n- dimensional 
subspace of T-L, and for 1 we let A„ the number defined as follows 

A„ sup inf (5.9) 

GeS"MGnV{F))s{0} {u,u)^ 
Then there exists a complete orthonormal system {ipn)n^ i such that 

Fipn = A„-0„, n ^ 1. 

In other words, the sequence {\n)n is the non- decreasing enumeration of the eigenvalues 
of F, each repeated a number of times equal to its multiplicity. Moreover, the sup in (j5.9p is 
attained for G equal to the span of {V'l, . . . , V'n}- 

Let us apply Proposition 15.31 to F = -A with domain 

P(-A) |n; n6C^(S), ^u(0)d0 = o|, (5.10) 

in H"^ (recall the definition of (H^^, (• , •)_|^) in Remark l2.2p and let us denote by £o{u,v) 
{u, -Av)_^ - 2tt fguv the Dirichlet form associated to -A. Note that Sq is well defined on 

D P(-A). Then, denoting by (aI"^^) n^i the sequence of eigenvalues associated to - A in 
Hi; 

^(-A) ■ r £0{U,U) 

XI, ' - sup mt ; — . 

Since the supremum is attained for G - {ipi, . . . , ipn} £ L^, one has in fact: 

sup 



Al-^).sup inf 



Ge5""e(GnL2)x{0} {u , u)_^ 

Secondly, note that one does not change the result by considering 1-Ai instead of -Ai. Hence, 
if ones denotes by £i{u,v) {u, (1 - Ai)v)_i g^^ the Dirichlet form associated to 1 - Ai, one 
deduces from [5l Eq.(2.47)] that £i is well defined on and that it is equivalent to Sq: there 
exists a constant C > such that 

Vu e L^, —£o{u,u) ^ £i{u,u) ^ C£q{u,u). 



Then, using again Proposition 15.3 



A^^^).sup inf ^^("'"^ 



sup 

Ge5"«e(G'nL2)\{0} [u , u}_^^^^^ 

Since the norms || • and || • ||-i,go are equivalent, one directly sees that there exist constants 
c,C > such that, for all n ^ 1 

cA(-^^^A(i-^^)^CA(-^). (5.11) 

One can prove similar bounds for A2 in the Hilbert space H^^ in the same way: ffist notice 
that any eigenvector which corresponds to a non-zero eigenvalue of A2 is necessarily with zero 
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mean-value, so that it suffices to work on the domain {v e L^, J^v - 0}. It is then easy to 
deduce from ()4.24p that both Dirichlet forms Sq and £2 (recall definition (I4.23P ) are equivalent 
on the subspace of with zero mean-value. In particular using Proposition 15.31 one easily 
obtains similar bounds as ()5.1ip for A2 and (15. 8p follows. □ 

Following the lines of [9], one can deduce that the linear operator Jq Tl(s)QTl{s)* ds is 
of class trace: indeed, it is easy to see from the definition (j4.2p of the perturbation B that B 
satisfies the assumption (5.59) in [9l p. 145], namely, i? is a continuous linear operator from 
into H~^, and there exists a constant c > such that for all /i e 2?, {Bh , h)-^^ ^ c || /i || jj^. 
Since Lemma 15.21 is true, the assumptions of [9l Prop. 5.25, p. 145] are satisfied, so that the 
operator /q Tl(s)QTl{s)* ds is of class trace. Then an application of |9, Th. 5.2] shows 
that the stochastic convolution Wl{-) is well defined as a predictable process in H^^. The 
assumptions of [9, Th. 5.4] concerning the existence and uniqueness a weak solution of (15. ip 
are satisfied and Proposition 15 . 1 1 is proved. 

5.3. Linear asymptotic behavior of the fluctuation process. We are in position to 
prove the main statement of Theorem 12.91 that is the asymptotic behavior ()2.19p of the mild 
solution ()5.7p . We will make a constant use of the spectral decomposition of the operator L 
(under the hypothesis of Theorems 12.51 and 12. 7p . We will refer in particular to the matrix 
representation (I2.17P of L. 

Note also that the continuous linear form idegi') can be represented, by Riesz 

representation theorem, as a scalar product w.r.t. some vector ( e H^^: 

4,,(-) = (C, •)h, • (5-12) 

The convergence ()2.19p is an easy consequence of Remark 12.81 and the following two proposi- 
tions: 

Proposition 5.4. The stochastic convolution WL^t) satisfies the following linear behavior, 
as t ^ +00: 0, where the convergence is in law. 

Proposition 5.5. For every initial condition rjQ, converges, as t ^ +00, to ip{r]o)deq. 

Let us now prove these two propositions: 

Proof of Proposition [5^ Recall that 1^ is a Q- Wiener process in Hq, which can be decom- 
posed into Ho = Span{dgq) ® G<o . Note also that the restriction on Hq of the projection 
Pq defined on H^^ by (12.16P coincides with idgq{-)dgq. With a small abuse of notations, we 
will use the same notation Pq for this restriction on Hq. Let us decompose the stochastic 
convolution into Wiit) - jQTL{t - s)Po dWs + j^Ti^t - s)P<q dWg, and treat the two terms 
separately: 

Let us consider the first term j jQTL{t - s)Pq dWs- One has successively, using that 
TL{u)dQq - doq for all u > 

^ f'TL{t-s)PodWs^^ rnit-sydogdWs, (5.13) 
u -J 1/ J 

- ^^do.Wt ~ (C, Wt)^^ , (recall (5.14) 
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Thanks to the Q- Wiener structure of W (see § 15. ip , one has 

(2 \ II II ^ II II ^ 

II i f^TUt - s)Po dW. \\J - (I (C , T^,)h, - {QC , Oh. , 

(5.15) 

which converges to as t ^ +oo. 

As far as the second term j^TL{t - s)E:o dWs is concerned, it is easy to see that it is the 
unique weak solution in Hq of 



= j\w,ds + P,oWt. (5.16) 



-t 

Wt 



Let us decompose evolution (j5.16p along this decomposition Hq = Span(5eg) © G<o : writing 

Wt ^ PoWt+ P<s) Wt := yt + zt, one has: 

I - f^PoLysds, 

[ yt ^ lo P<o LP<o Vs ds + P<oWt. 

Since the operator i^o LP^q has its spectrum in the negative part of the complex plane with a 
strictly positive spectral gap \k{L) and generates an semigroup of operators, it is immediate 
to see from the covariance estimates of stochastic convolutions (see [9l Th. 5.2, p. 119]) that 
there exist some to > and a constant c > such that for all t ^ Iq 

E(||y*|lHj' ^E(||yi||^J ^ce-^* (5.18) 

Consequently, one has 

E(NdlHj^ ^^(WPoLysW^^) ds^WdeqW^^ £B{\£9,giLy,)\) ds, 
^ b[\{C, Lys)iiJ\\deq\\-ti^) ds, (recall ^M) 
^ \\d9q\\n,\\L*C\\n,fJ^(\\ys\\n,)^^^ 

= II deq IIh, II L*C IIh, E (II IIhJ ds + £ E (|| ||hJ d.) . (5.19) 

In particular, it is immediate from (I5.18P and (I5.19[) to see that the following convergence 
holds: 

EflUtlln ) 

^ ^ -t^+oo 0. (5.20) 
Putting together (j5.18p and (j5.20p . Proposition 15.41 is proved. □ 

Proof of Proposition ] 5.^ Let us fix an initial condition 770 e H^^. Then X{t) ■- Ti^{t)riQ is the 
unique solution in of 



X{t)^r]o + Jj LXsds. (5.21) 



Decompose X{t) along the direct sum Go © G<o , that is X{t) = a{t)dgq + /3{t)p + Y{t), with 
Y(t) e G^Q. Then, projecting on d^q, p and G<o respectively (see (|2.17p ). one obtains that 
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(|5.2ip is equivalent to 

' - ^degiVo) + fo{m+ide<iiLP<oY{s))) ds, 

Vt>0, Pit) = tpivo), (5.22) 

, Yit) = P,oVo+IoP<oLP<oYis)ds. 

Then, since T/3g2.iJo(0 ^ semigroup of contraction whose infinitesimal generator has a 
strictly positive spectral gap XxiL), there exists in particular a constant c > such that 

Y{t) = Tp^oLP<oit)P<oVo and ||l"(t)||H^ ^ ce"^^* (in particular, j\\Y{t)\\^^ ^t^+oo 0). 
Then, using again (j5.12p . 

ait) _ id.qirjo) ^ , ^ , 1 r* 



t 



p{m) + ^ £^ae,iLP,oYis))ds, (5.23) 
iVo) + ^ (P<o ^"^C , >^(s))h, ds- (5-24) 



t 

Using the previous exponential bound for Y{s), it is easy to see that converges to ipirjo) 
as t ^ +00. The result of Proposition 15.51 follows. □ 

Appendix A. Gelfand-triple construction 

The construction of the weighted-Sobolev spaces defined in § 12.21 and used throughout the 
paper is based on the usual notion of Gelfand-triple (or rigged-Hilbert spaces) that we precise 
here. We refer to [H § 2.2] or [7| p. 81] for precise definitions. 

Namely, one can understand the closure of {/leC^(S); /§ ^ = O} w.r.t. the norm 
defined in ()2.7p as the dual space V' of the space V closure of {h e C\S); Jsh^O} with 
respect to the norm 

1 

Mv--{f^h'{e)'k{e)ddy . 

The pivot space is the usual L^(A), endowed with the Hilbert norm 

\h\\ ' 



|L2 



In particular, one easily sees that the inclusion V £ L^(A) is dense. Consequently, one can 
define T : L^(A) V by setting Th{v) - /g h{6)v{9) dO. One can prove that T continuously 
injects L^(A) into V' and that T(L^(A)) is dense into V' so that one can identify h€L^{X) 
with Th e V. Then for h e L^{X), 



\\h\\v, . \\Th\\v, . supif^ . \/T¥' (A-1) 

veV \\v\\v \ JS k 

where we used in (jA.ip Cauchy-Schwarz inequality for the lower bound and chose v' for 
the upper bound. This enables to identify H^^ with V'. 
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